Draft version March 13, 2013 

Preprint typeset using I^T^X style cmulatcapj v. 5/2/11 



BSG ALIGNMENT OF SDSS GALAXY GROUPS 

Zhigang Li 1,2 '*, Yougang Wang 1,1 , Xiaohu Yang 3 ' 4 ' 1 ", Xuelei Chen 1,54 , Lizhi Xie 2 ' 6,§ , Xin Wang 7 ^ 

Draft version March 13, 2013 

ABSTRACT 

We study the alignment signal between the distribution of brightest satellite galaxies (BSGs) and 
the major axis of their host groups using SDSS group catalog constructed by Yang et al. (2007). After 
correcting for the effect of group ellipticity, a statistically significant (~ 5a) major-axis alignment is 
detected and the alignment angle is found to be 43°. ± 0.4. More massive and richer groups show 
stronger BSG alignment. The BSG alignment around blue BCGs is slightly stronger than that around 
red BCGs. And red BSGs have much stronger major-axis alignment than blue BSGs. Unlike BSGs, 
other satellites do not show very significant alignment with group major axis. We further explore the 
BSG alignment in semi-analytic model (SAM) constructed by Guo ct al. (2011). We found general 
good agreement with observations: BSGs in SAM show strong major-axis alignment which depends 
on group mass and richness in the same way as observations; and none of other satellites exhibit 
prominent alignment. However, discrepancy also exists in that the SAM shows opposite BSG color 
dependence, which is most probably induced by the missing of large scale environment ingredient in 
SAM. The combination of two popular scenarios can explain the detected BSG alignment. The 
first one: satellites merged into the group preferentially along the surrounding filaments, which is 
strongly aligned with the major axis of the group. The second one: BSGs enter their host group more 
recently than other satellites, then will preserve more information about the assembling history and 
so the major-axis alignment. In SAM, we found positive evidence for the second scenario by the fact 
that BSGs merged into groups statistically more recently than other satellites. We also found that 
most of BSGs (80%) were BCGs before they merged into groups and earlier merged BSGs tend to be 
closer to their BCGs than other BSGs. On the other hand, although is opposite in SAM, the BSG 
color dependence in observation might indicate the first scenario as well. 
Subject headings: BSG - alignment - group of galaxies - SAM - SDSS 



1. INTRODUCTION 

In the current understanding of the hierarchical struc- 
ture formation scenario of cold dark matter (CDM), 
small dark matter halos form first from the anisotropic 
collapse of overdensities in the mass distribution, while 
larger halos grow by accreting surrounding material 
and/or by merging with other halos. It is well be- 
lieved that satellites in clusters or groups trace the dis- 
tribution of dark matter halos. Therefore, it is popular 
to study the shape or mass distribution of dark mat- 
ter halos using the distribution or kinematics of satel- 
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alignment signals between the distribution of satellites 
and the reference directions, e.g. the major axes of cen- 
tral galaxies, groups or halos, etc. 

Intrinsic alignment of galaxies relative to their 
host dark matter halos encodes plenty of infor- 
mation of galaxy formation and evolution his- 
tory and their su r round i ng filamentary structures 
(lAragon-Calvo et all l2007t iBailin fe Steinmetzl 120051: 



Lee et all 120081 ISlosar fc White! 12001 iVera-Ciro et all 
20111: iLibeskind et all 120121) . Observationally, three 
types of alignment have been extensively studied 
in the literatures: (1) the satellite alignment, be- 
tween the distribution of all satellite gal axies and 
the m ajor axis o f their central galaxy (iHolmbergj 
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20111): (2) the radial alignment, between the orientation 



of satellite galaxies and the central-satellite connection 
line or t he orientation of bri g htest central galaxies 
(BCG s ) (IPereira fc Kuhnl 120051: lAgustsson fc Brainerd! 
I2006bt iFaltenbacher et all 120071: lHao et alj 120111) : (3) 
the alignment between the shape of groups and the large 
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scale structure (iPaz ef^Lll2008l [201ll: iFaltenbacher et all 
I2009t IWang et all l2009f ). It is notable that the mea- 
surements of alignment between satellites and their 
host galaxies have a long and confused history. Some 
studies reported a minor axis alignment, while most 
of others preferred a major axis alignment . We refer 
those who are intereste d in th is history to lYang et all 
(|2006T) and IWang et all (l2008h for summaries. Now it 
is clear that the satellites preferentially lie along the 
major axis of host galaxies, i.e. a major axis alignment. 
Moreover, the alignment depends on the properties 
of the host halos, such as halo mass, color of central 
galaxies and satellites. For example, more massive 
halos usually show stronger alignment, and red BCGs 
and satellites show stronger alignment than the blue 
ones. These are consistent with the simulation results, 
where the distribution of subhalos is aligned with major 
axis of their host halos and mo re massive halo s show 
a stronger alignment signal dLibeskind et all 
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Faltenba cher et al. 1 1200a IKnebe et al.ll2008allbl . I201(£ 

Many dynamical processes can contribute to the align- 
ment of satellites associated with groups, which can be 
roughly divided into two classes - larg e scale environ- 
ment and the imp act of group potential (|Hartwickll200Ct 
IWang et al.ll20ldf ) . The filamentary structures surround- 
ing the group may readjust the satellites pointing toward 
the group. Additionally, the anisotropic accretion along 
these filaments will make the satellites preferentially dis- 
tributed along the group major axis. While, phase- 
mixing effect and relaxation will smear the memory of 
assembling history of satellites when they merged into 
the potential well of the groups. In this respect, more 
recently accreted or merged satellites contain more in- 
formation about the assembling history and should show 
stronger alignment with the major axis of the host halo 
than other satellites. 

It is expected that the brightest satellite galaxies 
(BSGs) enter the host group mainly through major merg- 
ers and play a determinant role in shaping the host group. 
If the BSGs enter the host halos statistically more re- 
cently, they will preserve the most information about the 
merger history. In this case, BSGs will be good tracers 
of the structure formation of their host groups. 

In this paper, we investigate the alignment of BSGs 
with the major axis of their host groups usin g the SPSS 
DR4 galaxy group catalog constructed by lYang et all 
(|2007| ). Different from the previous satellite alignment 
measurements, we focus on the alignment signal between 
the distribution of BSGs and the major axis of their host 
halos, not the major axis of the central galaxies. While 
the major axis of the central galaxy itself is misaligned 
with the project ed major axis of the host, e.g. at the 
level of - 23° (|Wang et all 120081 ) . Our observational 
measurements are then compared with the BSG align- 
ment in semi-analytic model (SAM) of g alaxy formation 
catalog constructed bv lGuo et al.l (|201lD . We also inves- 
tigate the distribution of the time when BSGs merged 
and compare them with those for other satellites in SAM. 

This paper is organized as following. In §2, we in- 
troduce the SDSS group catalog and SAM model briefly. 
The alignment estimator we used is described in §3. Then 
we show the results of BSG alignment in SDSS groups 
in §4 and in SAM in §5. In 55 we also show the distribu- 



tion of BSG merger time. Finally a brief summary and 
discussions are given in §6. 

2. DATA 
2.1. Groups of galaxies 

Our analysis is based on the SPSS DR4 galaxy group 
catalog, which is constructed bv lYang et al.l (|2~007l ) (here- 
after Y07) using an adaptive halo-based group finder 
(jYang et al.l 120051 ) . For a full description of the cata- 
log the reader is referred to their paper, but the major 
part of the catalog is presented briefly here. 

The base galaxy catalog is the New York University 
Value- Added Galaxy Cata logue of SPSS-PR4 (NYU- 
VAGC: lBla"nton et al.l[2005T ). As described in Y07, three 
group samples were constructed: Sample I, which only 
uses the 362,356 galaxies with measured r-band magni- 
tudes and redshifts from the SPSS, Sample II which also 
includes 7,091 galaxies with SPSS r-band magnitudes 
and redshifts taken from other surveys, and Sample III 
which includes an additional 38672 galaxies that fail to 
get redshifts due to fiber collisions and are assigned th e 
redshifts of their nearest neig hbors (jZehavi et al.|[2005h . 
Our analysis is based on Sample III. 

All magnitudes are e xtinction corrected (Schlcge l et all 
119981 ) and k-corrected (jBlanton et al.ll2003aF and evolved 
to rest- frame magni tudes at z = 0.1 usin g the evolving 
luminosity model of Blanto n et alJ ()2003bl ) . Stellar mass 
of galaxies ar e estim ated according to the fitting formula 
of iBell et all ((2003T) . The galaxies are assigned red or 
blue color accor ding to their bi-normal distribution in the 
01 (g-r) color (jBaldrv et alJ 120041 : iLi et al.l 120061 ). The 
galaxy is red if 01 {g - r) > 1.022 - 0.0651x - 0.00311a; 2 , 
where x = 01 M r — 5 log fa + 23. and blue otherwise 
(|Yang. Mo. k, van den Bosch! [20081) . The dark matter 
halo mass of group is estimated based on the ranking 
of the characteristic group stellar mass. Survey edge ef- 
fects have been taken into account by removing groups 
(~ 1.6%) that suffer severely from edge effects. Since 
the redshifts assigned according to the nearest neighbors 
are not very reliable, we discard groups whose brightest 
or second brightest galaxies have assigned redshifts. The 
total number of groups in the resulted catalog is 8,513. 

The brightest galaxy in the group is taken as the cen- 
tral galaxy and all others are satellite galaxies. In ad- 
dition, we also considered the most massive galaxy (in 
terms of stellar mass) as the central galaxy (MCG). As 
we have checked, the difference of the alignment signals 
between the luminosity (e.g. BCG) and stellar mass (e.g. 
MCG) indicators is too small to be noted. 

2.2. Semi- analytic model 

The semi-analytic catalog of galaxy populatio n used 
in this paper was constructed bv lGuo et al.l (|2011l ) based 
on the combination of Millennium (MS) and MS-II sim- 
ulation. The MS an d MS-II simulation are ru n by 
the Virgo Consortium (jBovlan-Kolchin et al. 2009]) and 
both contain 2160 3 particles. They adopted the same 
ACPM cosmology with parameters 57 t ot = 1, = 0.25, 
O b = 0.045, 0\ = 0.75, h = 0.73, a 8 = 0.9 and n s = 1. 
MS has a box size of 500 /i _1 Mpc and each particle has 
mass 8.6 x 10 8 /i _1 M , while the box size of MS-II is 
100 ft^Mpc and particle mass is 6.9 x 1O 6 /t. _1 M . Park 
matter halos in MS and MS-II simulation arc found by a 
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friend-of-friend (FOF) method with linking l ength being 
0.2 ti mes of the mean separation of particles (iDavis et all 
[1981 . The SUBFIND algorithm (jSpringel et all 120011) 
was applied to each group to identify all bound substruc- 
tures (subhalos). 

The galaxy formation models are implemented by al- 
lowing the galaxies to grow at the potential minima of 
the halos and subhalos in simulation. Each FOF group 
contains a central galaxy at the potential minimum of its 
main subhalo. Other satellite galaxies may reside at the 
potential minima of satellite subhalos, or be "orphans" 
which no longer correspond to any subhalos. Various 
semi-analytic effects arc considered, such as tidal torque, 
dynamical friction, star formation, gastrophysics, super- 
nova and active galactic nucleus feedback and galaxy 
mergers and etc. We refer those who are interested to 
IGuo et all (poTTh for more details about the SAM cat- 
alog. Here galaxies are also assigned red and blue col- 
ors according to the (g — r) color bi-normal distributions 
within these SAM galaxies. 

3. METHOD 

For a given galaxy group, the principal axis and its 
orientation on the sky can be calculated from the inertia 
tensor 



X, 



ATmom 

E 



Jr 2 ; i,j = 1,2 



(1) 



where N mom is the number of group members, (xj, n , Xj <n ) 
are the projected coordinates (with BCG as origin) of the 
nth satellite galaxy, and r is the distance between BCG 
and group members. The 1/r 2 weight is used to avoid 
overweighting on the group members at outskirts. The 
semi-major and semi-minor axes of the ellipse, L a and 
Lb, can be derived by solving the eigenvalue problem, 



An — L 2 X12 
A12 A22 — L 2 







(2) 



The direction of the major axis is given by the eigenvector 
r = [1, (L 2 — An)/Ai2], while the ellipticity, e, and the 
axis ratio, 77, are 



e = 1 - L b /L a , and 77 = L b /L a 



(3) 



We quantify the distribution of BSGs relative to the 
group major axes by the distribution function, 



P(9) = 



N(0) 

Ntnt 



(4) 



where 9 is the angle between the group major axis and 
the BSG-BCG connection line. The angle 9 is con- 
strained in the range 0° < 9 < 90°, where 9 = 0°(90°) 
indicates that the BSGs lie along the major (minor) axis 
of the host group. 

The strength of alignment can be quantified by the 
average of 9, 



(0) 



P(9)9d9 



(5) 



Then (9) < 45° indicates the alignment of BSGs with 
the group major axes, otherwise minor axes. The align- 
ment is stronger if (9) is further from 45° and there 
is no alignment at all if (9) = 45°. This estimator 



works well only if the groups arc spherical For elliptical 
groups, the alignment angle defined as this will generally 
depart from 45°, even if the BSGs are randomly dis- 
tributed. To consider the effect of group shape, one can 
generate mock samples in which the BSGs are randomly 
distributed according to the projected two dimensional 
axes of groups. To get the true alignment signals of the 
BSGs with respect to other satellite galaxies in the non- 
spherical groups, we define Pa(9) = P{&) — ^mock(^) + w 
with the normalization factor w, which can be calcu- 
lated by / P A (9)d9 = 1, and the alignment angle 9a = 
J P A (9)9d9 = (9) da ta ~(9) mock +45°. In practice, we con- 
struct mock catalog by randomly selecting the satellites 
as the BSGs in the groups, which reflect the distribution 
of all satellites. The ellipticity effect and shot noise are 
included automatically. The distribution functions for 
these mock samples are measured and averaged to get 
the mean angle, (#)mock, which ar e merely contributed 
by the shape of groups. We generated 1000 mock sam- 
ples to calculate the errors of distribution function and 
alignment angle. 

4. ALIGNMENT OF BSGS IN SDSS GROUPS 

In this section, we study the alignment signal between 
BSG-BCG connection line and the major axis of the host 
group. Firstly we measure the distribution function and 
alignment signal for the full sample in which each group 
has at least 6 members. The results are shown in FIG. 
[TJ The true alignment distribution function, Pa(9), is 
shown as black solid line with error bars come from 1000 
mock samples. A significant group major axis alignment 
of BSGs is detected. This can be seen from the shape 
of the alignment distribution function, and the resulting 
alignment angle which is 9a = 43°. ± 0.4 and smaller 
than 45°. 

As a comparison, the distribution functions, P{9), for 
data and mock sample arc also shown in FIG. [1] as green 
dashed line and red dot-dashed line respectively. Al- 
though the distribution function for data is much steeper 
and the corresponding mean angle, (#)data = 26°. 3, is 
much smaller, when the ellipticity effect is subtracted 
using mock sample, the true alignment signal is recon- 
structed. 

In the following we explore the dependence of BSG 
alignment on various group properties, including halo 
mass of groups in §4.1[ BCG and BSG color in §4 21 and 
group richness in §4.31 To answer the question whether 
the BSGs arc different from other satellites, we compare 
the alignment signal of BSG with that of other satellites 



4.1. Dependence on group mass 

Previous alignment studies in groups, which measure 
the alignment between the distributions of satellites and 
the major axes of central galaxies, have shown clear mass 
dependence, i.e., the alignment is stronger in mo re mas- 
sive groups (|Yang et all 120061 : iWang et all 120081 and ref- 
erences there in). In this section, we explore whether 
this mass dependence exists for BSG alignment. Limited 
by the size of SDSS group catalog, we divide them in to 
2 subsamples according to their halo mass. Low mass 
sample contains groups with mass log 10 Mh < 13.6 and 
high mass sample contains others. Where Mh is halo 
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Fig. 1. — Alignment of BSG and group major axis. The green 
dashed line shows the normalized distribution function of the angle 
between BSG-BCG connection line and group major axis. While 
the mean distribution function of mock samples is shown as red 
dot-dashed line. The black solid line with error bars is the true 
alignment distribution function. 
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Fig. 2. — Dependence of BSG alignment on group mass. Upper 
panel shows the results for low mass groups, and lower panel for 
high mass ones. 

mass with unit /i _1 Mq. 
are shown in FIG. [2j 



The BSG alignment for them 



4.2. Dependence on BCG and BSG color 

In this section we study how the BSG alignment 
depends on the colo r of BCG and B SG. Following 
lYang et all (|2006t ) and lWang etal] (120081 ). we divide the 



groups into four subsamples according to the color of 
BCGs and BSGs and measure the alignment signal of 
BSGs for them. The results are shown in FIG. [3] BSGs 
in groups with cither red BCGs or blue BCGs show sig- 
nificant alignment with major axis of their host groups, 
with a little bit stronger alignment in groups with blue 
BCGs. On the other hand, the results do show quite 
prominent BSG color dependence in that red BSGs show 
much stronger alignment with major axis of their host 
groups. 

4.3. Effect of group richness 

To study the effect of group richness, we construct sub- 
samples for groups with different richness and compare 
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Fig. 3. — BSG alignment dependence on BCG and BSG color. 
The upper panels show BSG alignment for groups with red (left) 
and blue (right) BCGs, while the lower panels for groups with red 
(left) and blue (right) BSGs respectively. 
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Fig. 4. — BSG alignment dependence of group richness. The 
results for groups with at least 6 members are shown in upper-left 
panel, 8 members in upper-right, 10 members in lower-left and 12 
members in lower-right panel respectively. 

the alignment signal for them. The results for groups 
with at least 6, 8, 10 and 12 members are shown in FIG. 
[H Although the measurement errors for rich groups are 
relatively larger, the trend that BSGs in richer groups 
are more strongly aligned with the major axis of their 
host groups, is clearly shown. 

4.4. Comparison with other satellites 

Finally, we study the alignment of other satellites with 
respect to the group major axis to see whether the BSG 
is special in this respect. In FIG. [5j we show the results 
for the second to fifth brightest satellites. Interestingly, 
all of them do not show any significant alignment, with 
alignment angles consistent with 45° at 2-a level. 

5. BSG ALIGNMENT AND MERGER TIME IN SAM 

In this section, we study the BSG alignment with group 
m ajor axis in SAM galaxy formation model constructed 
bv (Guo et all (1201 ID . To explore the origin of the BSG 
alignment, we also measure the time when BSGs merged 
into their host halos. 
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Fig. 5. — Alignment of other satellites for groups with more than 
6 members. Upper-left panel shows results for the 2nd BSG, upper- 
right and lower panels for 3rd, 4th (left) and 5th (right) brightest 
satellites respectively. 



5.1. BSG Alignment 

In SAM galaxy formation model, we select groups 
which have at least 6 members with stellar mass greater 
than 10 8 /i _1 Mq. The group major axis is estimated 
using the method discussed in Sj3] BSGs are selected 
by their rank of stellar mass in each group. The BSG 
alignment with group major axis in SAM and their de- 
pendence on group halo mass are shown in FIG. [51 
Clearly, BSGs show major-axis alignment and the align- 
ment angle is 42.7° ± 0.1. Furthermore, more mas- 
sive groups show stronger BSG alignment. Especially 
in massive groups with log Mh > 14.0, the alignment 
angle 37.7° ± 0.6 is much smaller than that obtained 
in SDSS observation. This can be understood in sev- 
eral ways. The fi rst possibility is that as discussed in 
iKang et all (I2007D . the various observational selection ef- 
fects the group catalog suffered from, such as interlop- 
ers (groups members do not actually belong to the same 
group) and incompleteness, tend to blur the true align- 
ment signal. Additionally, the hydrodynamical simula- 
tions, including baryon physics, have indicated that the 
condensation of baryons tend to make halo more spheri- 
cal or axisymmetric and cause the misali gnment between 
observations and N-bodv simulations (|Lin et al.l 120061: 
Debattista et alJl2008HPedrosa et al.ll2010tlTissera et al. | 
20101 : iBrvan et aJJ 12012b IBrvan et al.l 12013 iZemp et al.l 



Fig. 6. — BSG alignment with group major axis and its depen- 
dence on group halo mass in SAM galaxy formation model. 



| and references therein) . 



In FIG. we show the color dependence of the align- 
ment signals. First, for BCGs of different colors, the BSG 
alignments arc quite similar to observations and no sig- 
nificant color dependence is found. Second, for BSGs of 
different colors, quite similar to the observations, there 
exhibit prominent color dependence, however, in opposite 
direction. This discrepancy is quite interesting. As the 
stronger alignment of the BSGs with respect to the group 
major axes indicates the more recent accretion along the 
filaments, both results in SDSS observation and SAM 
are rather expected. In observation, the real Universe, 
these BSGs before their accretion into groups, arc pref- 
erentially located in filaments where the gas temperature 
might be quite high due to shock heating, etc., and are 
thus quenched and red. On the other hand, in SAM, al- 
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Fig. 7. — Dependence of BSG alignment on color of BCGs and 
BSGs in SAM. 

though these BSGs before accretion might also be pref- 
erentially located in the filaments, since SAM only make 
use of the halo merger trees where such kind of filament 
environment is not taken into account, and thus these 
BSGs (BCGs before accretion) might still be star form- 
ing and blue. 

Finally we show in FIGs. [5] and [HI the richness depen- 
dence of the BSGs and the alignment signals of other 
satellite galaxies. All these results are quite similar to 
those obtained from the SDSS observations, except that 
the second brightest satellite galaxies in SAM still show 
some weak major-axis alignment. 

5.2. BSG Merger Time 

Now we explore the time when BSGs merged into the 
host halos. We select all FOF halos with mass larger 
than 10 13 /i" 1 M© at z = 0. We look back the progeni- 
tors of BCG and BSG in each FOF halo along the merger 
tree, until we find the snapshot in which the progenitors 
of BCG and BSG belong to different halos for the first 
time. Then we define this snapshot as the BSG merger 
snapshot. We also get the merger time for other satel- 
lites (the second, third and fourth brightest in r band 
magnitude) using the same method. The distributions of 
merger time for the BSGs and other satellites are shown 
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Fig. 8. — Dependence of BSG alignment on group richness in 
SAM. 
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Fig. 9. — Alignment for other satellites in SAM. 

in FIG. [TOl Here one can see that there are more BSGs 
merged at low redshifts than other satellites. The aver- 
age time when BSGs merged is 9.87 Gyrs. While for the 
second, third and fourth BSGs, the average merger time 
are 9.58, 9.37 and 9.23 Gyrs respectively. These support 
the conjecture that BSGs merged into the groups statis- 
tically more recently than other satellites, although the 
signal is weak. Furthermore, 80.7% BSGs in their pro- 
genitor halos are BCGs before they merged with other 
halos. 

Finally, we check the scenario that earlier merged 
BSGs are closer to BCGs at redshift z = 0, by divid- 
ing the SAM catalog into 6 subsamples according to the 
time when BSGs merged and calculating the distribution 
function of distances between BCGs and BSGs at z = 0. 
As shown in FIG. [Til most of earlier merged BSGs are 
located within 0.2 Mpc to their BCGs. While only a tiny 
fraction of later merged BSGs are located within 0.2 Mpc 
to their BCGs and most of them are far from the BCGs. 

6. SUMMARY AND DISCUSSIONS 

In this work, we studied the alignment of BSG-BCG 
connection line with group major axis using the SDSS 
DR4 group ca talog const r ucted by Y07 and the SAM 
constructed bv lGuo et al.l (|2011l ). We found good agree- 
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Fig. 10. — The distribution of merger time of BSGs and other 
satellites. The black solid line shows the distribution of number of 
BSGs which merged into groups with respect to redshift or time of 
universe. The magenta dotted, green dashed and red dot-dashed 
lines are results for the second, third and fourth brightest satellite 
galaxies, respectively. The average time when mergers happened 
are 9.87, 9.58, 9.37 and 9.23 Gyrs for the BSGs, second, third and 
fourth BSGs respectively. 
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Fig. 11. — Distribution of distances between BCG and BSG at 
z = in subsamples of groups whose BSGs merged in different 
time. The black solid line are for groups whose BSGs merged at 
the time between 2.16 to 3.52 Gyrs, while magenta dotted line for 
BSG merger time between 3.53 to 5.33 Gyrs, and so on. 

ment of major-axis alignment of BSGs in observations 
and SAM. We also studied the distribution of merger 
time for BSGs and other satellites in SAM and found pos- 
itive evidence to support the scenario that BSGs merged 
into their host groups more recently than other satellites, 
and preferentially along the filaments. The main results 
are summarized as following: 

• A significant group major-axis alignment of BSG 
is detected. The alignment angle is 43°. ± 0.4. 
The confidence level of BSG alignment detection 
reaches 5er, owing to the large group catalog. BSGs 
in more massive groups and richer groups show 
stronger alignment with the major axis of their host 
groups. Furthermore, the BSG alignment around 
blue BCGs is slightly stronger than that around red 
BCGs. And red BSGs have much stronger major- 
axis alignment than blue BSGs. 
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It is worth to distinguish our BSG alignment from 
other satellite alignment measurements. Most of 
them reflect the global distribution of all satellites 
along the major or minor axes of groups, while we 
measure the preference of BSG distribution along 
the major or minor axes relative to the other satel- 
lites. 



• Satellites other than BSGs do not show any signif- 
icant alignment. This may indicate that the BSGs 
are quite different from other satellites in the for- 
mation history of groups. 



• In SAM, BSGs are found to have major-axis align- 
ment with similar strength as in observations. The 
alignment show strong dependence on group mass 
and richness in the same way as observations, ex- 
cept that this dependence in SAM is stronger. 



• Discrepancy also exists in that the SAM shows op- 
posite BSG color dependence, which is most proba- 
bly induced by the fact that SAM does not contain 
any large scale environment ingredient. 



• BSGs in SAM are found to be merged into the host 
halos statistically more recently. And most of the 
BSGs (80%) were BCGs before their merger and 
earlier merged BSGs are closer to their BCGs at 
redshift z = 0. 
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The detected alignment of BSGs along the group major 
axes and no alignment for other satellites can be under- 
stood if the two popular scenarios of galaxy formation 
are both true. The first scenario is that the accretion of 
satellites are preferentially along the filamentary struc- 
tures surrounding the dark matter halos of host galax- 
ies. The second one is that most massive satellites are 
accreted into the host halos more recently and so they 
have experienced less phase-mixing and relaxation and 
better preserved the memory of their accretion history. 
Using SAM galaxy catalog, we find positive evidence to 
support the second scenario. To fully understand the 
detected alignment of BSGs, detailed analysis of the ac- 
cretion history of BSGs in SAM are necessary, which is 
beyond the scope of this paper. However, although is op- 
posite in SAM, the BSG color dependence in observation 
might indicate the first scenario as well. 
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